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1. INTRODUCTION 

The key tasks to be performed under contract No. NAS9-16873 were: 

(1) Given the reflectance data for a crop canopy and a canopy reflec- 
tance (CR) model which satisfactorily represents the crop (e.g. the Suits 
and the SAIL model for a uniform and homogenous crop canopy), determine 

how correctly and uniquely can one determine the canopy parameters of import- 
ance, and how sensitive is this determination to the variation in the reflec- 
tance data? 

(2) Implement on a computer the inversion of a CR model and use it to 
determine the agronomic parameters using measured data on the reflectance of 
representative canopies. A comparison of the calculated parameters with the 
field measured parameters will then provide an assessment of the possibility 
of estimating agronomic parameters from the CR data via inversion of a CR 
model . 

Work was done towards both the tasks. The reports on progress have been 
periodically presented to the NASA-managements and at various professional 
meetings. Since the detailed results have also been documented in the form of 
papers which have been published in the open literature, in this final technical 
report we will only highlight the activities carried out and the progress made 
and summarize results obtained to date. We refer the readers to various docu- 
ments listed in Appendix A for further details. (In Appendix B, we have added 
copies of the first pages of the reports/papers which have been published). 

2. HIGHLIGHTS OF THE ACTIVITIES 

The principal investigator and his associates focused their efforts in 
developing a technique for inverting CR models, i.e. a technique for estimating 
all the canopy parameters, occur'ng in a CR model, from the measured canopy 
reflectance data. In the published report Rl, v/e laid down the strategy for 
our efforts. According to it, we first chose a simple CR model - namely, the 
Suits' model for a one-layer, one-dimensional homogenous canopy. This model 


uses six canopy parameters; 
Agronomic Parameters: 
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Hd, Vd 


Spectral Parameters; P>T,Pg 

Illumination Source Parameter: SKYL 

iiere H and V, respectively, are the projections on the horizontal 
and vertical planes of the average leaf area per unit volume, d is the 
height of the canopy, and SKYL is the diffused fraction of the incident 
solar radiation, p is the leaf reflectance, t is the leaf transmittance 
and Pg is the soil reflectance. The dependence of reflectance on v/ave- 
length X is Implicit in the dependence on X of p,x,p , and to some extent 
of SKYL. 

In addition to the above canopy parameters, the CR also depends upon 
the illumination and view directions parameters (0^, 0^, \p). Here 0^ is 
the solar zenith angle, 0^ is the view zenith angle, and is the relative 
azimuth angle (0^ ^ ^ leo**) between sun and view directions. 

Initially for the purpose of investigating the invertibility of this 
simple Suits model, v/e used "perfect” or "error free" CR data. That is, 
we used the Suits', Model to calculate the CRs for a set of values of the 
canopy parameters. Then these reflectances v/ere treated as measured values 
and used in inverting the CR model to obtain canopy parameters. The values 
of these parameters were compared with those used in the original forward- 

direction calculation. 

« 

The basic procedure for inverting a CR model is as follows. Let R^ be 
the CR computed for a set of canopy parameters for the so'lar/view direction 
(0- (i}» 0„ {i)» and ij; (i). Let R^ be the measured CR for the same set of 

V f 

solar/view angles. We define a merit function F by 
n , 

F = E w. (R. - R'.r 
i=l ^ ^ T 


( 1 ) 
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Here, the sumnation 1s over the number of solar/view angles, and 
the canopy parameters appear as unknowns in R^. The W|'s are weight 
factors which could all be set equal to 1, or be given unequal values 
to reflect the relative importance and accuracy of the various obser- 
ved CRs. In the absence of any other information, we set all these 
weight factors equal to 1. The unknown canopy parameters which best 
match a given set of values, R^. ', of observed CRs, are those for which 
the function F is minimal. (Hopefully this minimum will exist and be 
unique.) 

The basic procedure for finding the (global) minimum of F is to start 
with an initial guess for the unknown canopy parameters, to calculate F 
using the CR model, and to use some method to choose successive values for 
the parameters until the computed F takes on its minimum value. 

Initially we used a general purpose nonlinear optimization procedure 
to find the minimum value of F. This procedure was designed to ensure that 
one will always reach the global minimum of F while iterating on the canopy 
parameters. Using this procedure, we showed that the agronomic parameters 
of the Suits' model Hd, Vd (which are related to the leaf area index by 
LAI e d(H'^ + V )'®), could be estimated from the CR data in one v/avelength band 
(preferably near-infrared band) for a few solar/view directions, provided that 
the o!,Iidr canopy parameters (p,T,pg and SKYL) are known (i.e. kept fixed at 
their '' lues in the optimization procedure). 

To assess if canopy parameters other than LAI could also be estimated from 
the CR data through inversion of the CR model, v/e also chose p and t as unknown 
parameters. With these four parameters as unknown, we Initially found that, 
depending upon the initial guess for the canopy parameters, one may or may not 
succeed in determining their values by the inversion process. That is, for 
certain initial guesses one obtains values which are correct, while for others 
one is unable to obtain correct values. Instead a so-called slow convergence 
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problem is encountered in which unrestricted computer time is required 
to succeed in obtaining the correct’ values. 

We analyzed the slow convergence phenomenon and developed a so called 
'angle transform' approach v/hich alleviated the problem somewhat. In this 
approach, one constructs combinations or transforms of CRs at various solar/ 
view angles which are either sensitive or insensitive to a given agrophysical 
variable. Ideally, one would like to have transforms which are minimally 
sensitive to those variables/parameters which are prone to error in measure- 
ment and maximally sensitive to those which one wishes to estimate by the 
inversion, technique. Then one could sv/itch from one transform to another to 
estimate one set and then another set of canopy parameters. For example, 
we were able to define an angle transform which was very sensitive to ptr 
and relatively insensitive to individual values of p or t and other transforms 
which are sensitive to H/V and p/t. Using these angle transforms we showed 
that the slow convergence problem could be avoided and we could determine all 
the four canopy parameters Hd, Vd, p, and t, from the CR data. The results 
of our efforts are described in detail in published report R3. 

In this report v/e also presented an analytical procedure for assessing 
the levels of errors in the estimation of canopy parameters as a function of 
errors in the measurement of CRs or, equivalently, as a function of accuracy 
v/ith which a CR model represents the canopy reflectance. This procedure 
shov/ed that if p and t are kept fixed in the inversion procedure, the impor- 
tant agronomic variable LAI could be estimated quite accurately. However, 

• . v/hen these tv/o parameters are treated as unknown, the estimation of LAI could 
be quite erroneous. Further the level of error depends upon the solar/view 
angles used in the inversion process. In fact, the error analysis could be 
used to define optimal set of solar/view angles. 

The angle transform approach, though addressed the slow convergence pro- 
blem, did not satisfactorily alleviate the problem; the convergence to the 
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optimum values for the canopy parameters still took several thousands of 
Iterations. 

Instead of using a general purpose optimizer, v/e then developed an 
optimizer especially for this problem which is much more efficient, re- 
quiring tens rather than thousands of iterations. This optimizer is des- 
cribed in detail in the published report R4. We shovved that, with this 
procedure, one can, in principle, estimate all six canopy parameters appear- 
ing in the Suits model for homogenous single layer canopies, using only CR 
data. In other words, we showed that such a Suits model is mathematically 
totally invertible. We also carried out an analysis of the accuracy of the 
estimation of canopy parameters as a function of random errors in the CR 
data using the techniques developed in report R3. We investigated how this 
accuracy depends upon the nominal values of the canopy parameters, the num- 
ber of unknown parameters, and tho number of solar/view angles for which the 
CRs are used in the inversion process. 

We then applied the techniques for CR model inversion and for analyzing 
the accuracy of estimation of canopy parameters to two more complex models - 
the SAIL (scattering by arbitarily inclined leaves) model due to W. Verhoef 
and N. Bunnick of the National Aerospace Laboratory of the Netherlands and the 
CUPID model due to J. Norman of the University of Nebraska. Suits model 
assumes that the leaves are either horizontally or vertically inclined while 
both of these models allow arbitrary leaf angle distributions (LAD). Both 
of these models use fraction of leaves at discrete leaf inclination angles 
as parameters. To minimi ■'e the number of unknown canopy parameters, we looked 
for a simple continuous distribution (v/ith minimal number of parameters) v/hich 
could represent the leaf angle distribution well. In published report R8, we 
showed that a simple’ beta distribution, characterized by two parameters, )i and 
V, represents well not only the v/eTl~known leaf angle distributions (planophile, 
plagophile, erectophile, extremophile, uniform and spherical) but also measured 
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distribiitions for soybean, wheat and blue grama grass canopies. With 
this simplification, the number of canopy parameters, for both the 
SAIL and the CUPID models, increased only to 7 from 6 for the Suits' 
model. We showed that both SAIL and CUPID models are also mathemati- 
cally totally invertible. The details of the inversion and the error 
analysis are given in published reports R2 and R5. 

The error analysis 'for-' aTl- these’ iriodels**shoVJed' that' when one es^timates 
the canopy parameters using only the bidirectional ’CR data, 'for certain 
nominal' values of the canopy paramatersj'the estimated value of a parameter 
can change by several thousand percent v/hen CRs are randomly ch.anged by 
1 %. This analysis led us to conclude that, in general, for the expected 
levels of errrors in the measurement of CRs and the accuracy with which 

* ■* I 

these models are likely to represent CR, one can not practically determine 
the agronomic parameters using CR data alone . Such a detemi nation , in 
general, will require ancillary data e.g. on the reflectance and transmi- 
ttance of vegetation elements and on the soil reflectance. 

The activities and the results summarized above were all directed 
towards completing task (!) mentioned in the Introduction. Towards complet- 
ing the task (2), we implemented a computer program for carrying out the in- 
version of a CR model. This program is general enough to be usable with any 
CR model. Using the SAIL model, we tested if the inversion technique can be 
used to estimate LAI and LAD using measured CR data. In published report 
R6 we showed that for a fully covered soybean canopy, one can indeed deter- 
mine fairly accurately the LAI and the average leaf inclination angle (ALA) 
using measured CR data (collected, using ground based sensors, for about 
50 view angles over a period of about 15 minutes) and easy to measure canopy 
parameters (leaf reflectance and transmittance, soil reflectance, and frac- 
tion of diffused skylight). The estimated LAD appears to be somewhat erro- 
neous; it is a narrov/er distribution (peaked at an angle close to the correct 
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value of ALA) than the observed one. 

To determine if the inversion technique can be used to estimate 
LAI for an inhomogenous canopy such as that for a forest using remotely 
sensed bidirectional CR data and some (ground) measured ancillary data on 
spectral and architectural parameters of the canopy, we applied the tech- 
nique for a black spruce canopy for which the CR data was collected using 
a C-130 aircraft, equipped with a Thematic Mapper (TH) sensor. In publi- 
shed report R7 we showed that in principle LAI can be estimated. However, 
because of the complexity and inhomogenity of the black spruce canopy, it 
can not be accurately determined with just 7 bidirectional CRs which v/ere 
available. Ideally, one would require measurements of CRs for more solar/ 
view angles. In addition, some improvements on the CR model (e.g. inclusion 
of shadow effects) so that it represents more accurately the CR of a forest 
canopy are needed. 

We also briefly carried out an analysis of optimal solar/view angles 
for measuring CR for v/hich the estimation of LAI is likely to be most ac- 
curate. The results of this analysis were presented in an oral report 
(PIO) when specific recommendations were made for CR data collection stra- 
tegy. This strategy was implemented for the field v/ork carried out during 
the summer of 1984. Whether the concept of optimal solar/view angles is 
valid or not must await the analysis of this data. 

3. FUTURE RECOMMENDED ACTIVITIES 

We successfully accomplished the two tasks specified in the contract. 

The research on the inversion of the canopy reflectance model, which we 
carried out is obviously not yet complete. We recommend that the following 
activities, divided into three major categories, all related to the inversion 
of the CR models be pursued. 

(1) Inversion of the One-Dimensional CR Models 

The activities in this category could be simply described as the natural 
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continuation of the work already done and summarized in the preceding 
section. They are essentially of two types: • 

(a) Additional apsects of Inversion and sensitivity analysis of the 
existing CR models investigated so far (Suits, SAIL, and CUPID models). 

(b) Modification of these models to make them more realistic and 
then the inversion and sensitivity analysis of modified models. 

In the first type, we recommend that the following activities be 
carried out: 

{ a . 1 ) Use of CR Data in Many Wavelengths or Spectral Bands 

As the number of wavelengths or spectral bands increases, the data 
base for estimating LAI and LAD increases. Thus it can be argued that one 
could substitute wavelengths for solar/view angles and obtain these agronomic 
parameters with fewer bidirectional CR measurements. The advantages and dis- 
advantages of using many wavelengths over many solar/view angles should be 
quantified. Here, it should be noted that as the number of wavelengths in- 
creases, the number of ancillary parameters (the spectral parameters like leaf 
reflectance and transmittance, soil reflectance and fraction of diffused sky- 
light) also increases. 

The approach for this quantification is straightforward. It v/ill involve 
increasing the number of canopy variables and the number of R'^ in Eq. (1), 
where now the index i will refer to solar/ view angles as well as wavelength/ 
band. By trading off CR for various angles with CR for various wavelengths, 
hope to arrive at an optimal combination for which LAI and LAD estimation 
Is most accurate. This should then be tested by using available CR data sets 
on various vegetative canopies. 

( 3 • 2 ) Use of Linear and Non-Linear Transfoniis of CRs for Various Solar/ 
View Angles and Various Spectral Bands 

In the present inversion technique, the measured CRs are fitted to those 
calculated from a model in a least square fashion. It is not obvious that the 
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choice of such a fit Is the best. In other words, If in Eq. (1) one 
replaces by a linear or non-linear function of several R.|S, for 
different solar/view angles and for different spectral bands, could 
the sensitivity of F to LAI and LAD change so as to make their esti- 
mations more accurate? The basis for raising such a possibility Is 
the well known linear {greenness/brightness) transform of CRs for 
four MSS bands or seven TM bands. Greenness is very sensitive to LAI 
and not to soil reflectance while brightness Is very sensitive to soil 
reflectance and not to LAI. Another simlllar reason Is the concept of 
angle transfonn mentioned In the preceding section. 

It is recommended that the Inversion and sensitivity analysis using 
various linear and non-linear transforms of CRs for various solar/view 
angles and for various Wavelength bandiS 'je carried out. In fact, sensi- 
tivity analysis could be used to search for such transforms. The goal 
should be to choose that form of F which has minimum sensitivity to an- 
cillary variables/parameters (which may have some errors In their measure- 
ments) and maximum sensitivity to those (LAI and LAD) which one wishes to 
estimate by the inversion technique. 

(a. 3) Inversion of Radiance Data Inside the Canopy 
The LAI estimation based upon CR measurements above the canopy Involves 
the use of a model which accurately represents the interaction of downwelling , 
radiation as well the interaction of upv/elling radiation, reflected from back- 
ground, with the canopy elements. On the other hand, If one could measure the 
radiation flux inside the canopy at various locations and relate It (through 
essentially "half" of the CR model) to LAI and LAD, it appears that the esti- 
mation of LAI and LAO could be simpler and more accurate. The technique for 
below canopy measurement of sunfleck area for estimating LAI, developed by 
John Norman, essentially makes use of this simplicity. 

It 1s reconinended that the fluxes inside the canopy at different heights, 
using the simple CR models discussed above, be calculated and then the Inver- 
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slon and sensitivity analysis with these fluxes as the measured variables 
rather than CRs be carried out. This analysis should then be used to 
quantify the accuracy of estimation of LAI and LAD for typical canopies 
and for specifying the most desirable set of measurements (Incident solar 
direction, location of sensors In the canopy, etc.). This work should 
lead to a rapid method for ground measurement of LAI and LAD which are 
heretofore labor Intensive to measure. 

(a. 4) Optimal Solar/View Angles for LAI and LAD Estimation 
On the basis of sensitivity analysis, we have conjectured that there 
are certain solar/view angles for which LAI and LAD estimation Is likely to 
be most accurate. The set of these angles varies somewhat with the nature 
of the canopy but still one could provide some guidelines for selecting these 
angles for all type of canopies. 

We recommend that one should look Into details of the Interception of 
radiation by the canopy elements, as Imbedded in the above mentioned CR models, 
as a function of solar/view angles, We expect that one will find that for the 
optimal angles the gradient of Interception as a function of LAI and LAD Is 
maximum. This enquiry, as a function of values of ancillary canopy parameters, 
should provide a fundamental basis for optimal selection of solar/view angles 
and wavelengths for the purpose of estimating agronomic variables from the 
*'R data. 

Parenthetically we note that If one could identify a small set of optimal 
^Olar/vlew angles for which all the canopy parameters can be determined rea- 
sonably accurately from the CR data, then these parameters can be used to 
calculate bidirectional reflectances for ^ solar/view angles. These reflec- 
tances can then be integerated to find the hemi -spherical reflectance and the 
vegetation albedo. 

(a. 5) Angle Transforms for Filtering out Terrain Slope Effects 
In the analysis of remotely sensed CR data It Is usually assumed that the 
topography of the scene 1s flat, l.e. , all surfaces or materials are assumed 
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to be viewed at same sensor and illumination angles. In mountainous 
terrain, however, the local surface normal varies and, consequently, 
a wide range of effective view and illumination angles are obtained, 
even for a single image. In these cases the interpretation of CR data 
could lead to erroneous agrophysical conclusions. 

Several authors have tried various correction procedures to improve 
the results obtained. These approaches arc either essentially empirical, 
making no assumptions concerning the physical behavior of scene elements, 
or utilize some assumed natural characteristics of the scene element. 

It is recommended that the angle transfonn approach be used to develop 
a procedure which is minimally sensitive to local surface normal. One 
specific approach is to characterize the normal to the ground surface with 
two angles {zenith and azimuth) and calculate the CR reflectances for various 
illumination/ view angles by using one of the above mentioned CR models. Such 
a calculation is straightforward and is essentially equivalent to using a 
different set of illumination and view angles. One can then seek transforms 
of CRs for- these angles which -are insensitive to the angles defining the normal 
to the surface. If such transforms can be found (preliminary calculations 
support such a possibility), then one would have a simple technique to filter 
out the effects of terrain slope. 

We now discuss the recommended activities of the other type involving 
modifications of the one-dimensional models constantly referred to in the 
preceding description of recommended future activities. All these activities 
are designed to make the CR models more realistic and thus to ensure that the 
estimated values of agrophysical parameters, based on inversion of these models, 
is more accurate. The following three activities are recommended for future 
investigation. 

(b.l) Inclusion of Specular Reflectance 

Current CR models do not adequately account for specular reflectance, 
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(dlrect mirror-tike reflection from first-surface leaves), though it has 
shown to be quite significant for a number of vegetations and it has been 
known for some time that leaf reflectance is quite directional. Since the 
model inversion technique relies on angular variation of the CR data, more 
accurate representation of the angle dependence is crucial for its success 
with actual field data. 

It is recommended that a model for specular reflectance, based on 
geometrical optics be developed. This approach is knov/n from the physics 
literature to work for reflectance from "rough" surfaces. In this approach 
one will consider leaves as reflecting "facets" of a surface, distributed 
according to the leaf angle distribution of the canopy. The angular depen- 
dence of the specularly reflected radiation will thus be a "trignometrical 
convolution" of the direction of the incident radiation and LAD. The total 
reflectance could then be decomposed into the sum of a direct specular part 
(determined by geometrical optics approach) and a uniformed diffused part 
(determined by using one of the CR model referred above). Since electro- 
magnetic waves combine according to the superposition principle, this is a 
very natural and well-justified decomposition. 

The modified model should then be subjected to the inversion technique 
v/ith measured CR data as the input, to determine if the new model does indeed 
give more accurate values of the agronomic pa>"dmeters. It should be pointed 
out that the relative proportions of specular and diffuse components will be 
taken as an unknown parameter in the inversion process. Its values can then 
be calculated with the independent measurements of such proportions carried 
oi't by Vanderbilt and Grant at LARS using polarized light, 

(b.2) Inclusion of Angular Dependence of Leaf Reflectance and Trans- 
mittance 

In the estimation of LAI and LAD, from the measured bidirectional CR data 
via Inversion of a CR model, the values of leaf reflectance p, and transmit- 
tance, T, used are usually measured with incident radiation located in the 


nadir direction. Such measurements should tend to give a lower value for 
p and a higher value for t than will be felt by the incident radiation in 
a, typical canopy. Moreover, this differential will be dependent upon the 
direction of solar radiation, leading to a modification of the dependence 
of the canopy reflectance on the i 11 unii nation direction. 

It is recommended that the angular dependence of p and t in the CR 
model be included and the modified model be investigated to determine if 
the estimated values of the agronomic parameters through the inversion tech- 
nique improve. The specific dependences to be tried could be guided by the 
direct experiments on the angular dependence of p and t on the incident and 
view angles and by various models of leaf reflectance and transmittances. 

(b.3) Inclusion of Shadowing Effects 

The CR, in the near infrared bands, for low density trees stands has 
been ."ound to be higher than for high density tree stands rather than the 
other way around, as implied by the simple canopy reflectance models used 
so fa” in the inversion technique. When these models are used for low den- 
sity tree stands, as to be expected, the estimated values of agronomic 
parameters are quite erroneously high (see published report R8). It is 
believed that the higher CR for low density tree stands is due to shadow 
effects (shadow of a tree on itself, on other trees and on ground). 

It is recommended that these simple models be modified to include the 
shadowing effects. An approach could be geometrical, with trees assumed to be 
spatially distributed according to a pre-assigned distribution. The amount 
and nature of shadowing will be dependent on the illumination direction and 
the effects of shadowing on the measured CR will be dependent on the view 
direction. 

We now describe the recommended future activities. 

(2) Inclusion of Atmospheric Scattering 

In all the inversions of canopy reflectance models we have neglected the 
effect'; of atmospheric scattering of reflected radiation. This simplification 
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vas not detrimental because we used CR only in the near infrared band. 

Also, the CR used in the inversion process were either measured on the 
ground and/or in very clear sky conditions. For remotely measured CRs, 
especially in the visible band, the atmospheric effects are obviously 
very important and for certain atmospheric conditions may in fact domi- 
nate. As mentioned above under the discussion of optimal solar/ view 
angles, when atmospheric scattering is neglected high view zenith (and 
solar zenith) angles are more desirable than low zenith angles, to maxi- 
mize the accuracy of estimation of agronomic parameters like LAI and LAD. 

For such angles, however, due to increased radiation path length between 
canopy and the sensor, the atmospheric scattering effects should tend to 
be maximum. Therefore, before one could develop an optimal strategy for 
bidirectional CR measurements one must include atmospheric scattering 
effects in the model to be inverted or develop methods to "filter out" 
these effects. 

There are two basic approaches to include the effects of atmospheric 
scattering. One is to treat atmosphere and the vegetation canopy as a single 
system and solve the radiative transfer equation for complex boundary con- 
ditions. Such an approach is being pursued by Gerstl and Diner and Martonchik. 
The other approach is to use an atmospheric scattering model which generates 
parameters which can be used in the CR models to couple their outputs to the 
atmospnere. (Parameters needed are atmospheric transmittance from nadir, 
path radiance, and the direct and diffuse components of solar flux incident 
on the earth's surface as a function of v/avelengths and solar zenith angle.) 
Such a model then needs to be incorporated into the vegetation CR model. 

The first type of approaches are obviously more rigorous, representative 
and elegant. However, they are still under development anu not be gene- 
rally available for a few years. Also, it appears that computet* time for 
solving the radiative transfer equation may be excessive enough to preclude 
it from being an attractive candidate for a CR model to be used in the inver- 
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sion process (Recall that inversion techniques uses an iterative procedure 
and in each iteration the CR model is used to calculate CRs for various 
solar/view angles, for a given set of canopy parameters). 

We therefore recommend that the second approach be pursued using 
atmospheric scattering model, such as due to Dave and/or the LOWTRANS 6 
model developed at Hanscom Air Force Base in Massachusetts. These models 
should then be used to calculate atmospheric scattering for a few model 
atmospheres that range from "clear” to "murky".. Aerosol distributions could 
be limited to continental spring, continental summer and maritime summer. 

For each of these distributions, one could choose say, three optical thick- 
ness, from low to high. 

The incorporation of the atmospheric model with the vegetation canopy 
reflectance models can be done in the following two v/ays; 

(1) Coupling an atmospheric scattering model with a CR model : One way 
this can be achieved is by extending a multi-layer canopy model by adding 
two or three layers atmospheric model at the top. One approach would des- 
cribe atmospheric parameters in terms of CR model parameters (e.g., scattering 
cross section expressed in terms of effective LAI and LAD). Another approach 
simply couples existing atmospheric scattering and CR models so that the lower 
boundary conditions for the atmospheric model match with the upper boundary 
conditions for the canopy model. It should be noted that since either approach 
leads to a combined model whose upper boundary condition is specified by the 
solar energy curve, a simplification in parameter specification is achieved 
namely, it will not be necessary to specify fraction of incident dir'^used 
radiation. 

(2) Treating atmospheric scattering and canopy reflectance models as 
uncoupled; Here the secondary interactions between various components 

of the total remote sensing system (canopy, atmosphere, sunlight) are essen- 
tially ignored and the approach is therefore considerably simpler. The 
atmosphere and the canopy are treated as separate entities which are coupled 
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only in that flux from the canopy passes through the atmosphere. The 
scattering of flux back from the atmosphere to the canopy is treated only 
in a simplistic way without considering detailed interactions. Likewise, 
the atmosphere contributes to the canopy only in that it scatters incident 
sunlight to provide diffuse radiance on the canopy. 

Both of these approaches should be explored. The CR model coupled 
with atmospheric scattering effects included should then be subjected to 
the inversion and sensitivity analysis in the same fashion as we have done 
so far. One could also investigate the possibility of using an angle trans- 
form to "filter out" the effects of atmospheric scattering. 

Finally, we now describe the recommended research activities in the 
third major category. 

(3) I nversion of Complex Canopy Reflectance Models 

Simple CR models represent fairly well only uniform vegetation canopies 
but are generally inaccurate for canopies with complex geometries such as 
for crops planted iri rows, and for trees. For such canopies, one must resort 
to complex models. 

It is recommended that the inversion and sensitivity analysis of a few 
such complex models be carried out. The choice of the models will be depen- 
dent upon the maturity of a model and the availability of the softv^are imple- 
menting the model. 

Suits row model is one such candidate. Not only this model is mature 
and a software package is available, it is simple enough that it can be 
hiodified and made more realistic in the same v/ay as described above for the 
other homogenous models. 

Norman and Wells three dimensional model is another potential candidate. 
Kimes and Kirchner model is still another possibility. (It is somewhat less 
developed than the Norman and Wells model). 

There are other models which are actively being developed. It is quite 
conceivable that one of these models could be chosen over those mentioned 
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above, depending upon the advancement of the model and the v/illingness of 
the author(s) to make the softv/are available to other investigators. 

The inversion technique, with these complex models, should then be 
tested using measured CRs of heterogenous canopies such as those for forests. 
Such data is being collected by investigators at many institutions (e.g. 
NASA-Goddard Space Flight Center, NASA-Johnson Space Center, Purdue Univer- 
sity - LARS, Oak Ridge National Labs.). 
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APPENDIX A 

REPORTING OF THE PROGRESS AND PUBLICATIONS 


In this appendix we v/ill list the reporting of the progress made 
throughout the tenure of this contract. The reporting was done both 
in the form of briefing as well as in the form of technical reports 
and papers by the principal investigator as well as by other personnel 
associated with the contract. The copies of the slides used in the 


briefing and reports/papers have already been provided to the Technical 
Monitor. It should be noted that some of the technical results were 


only partially supported by this contract. 


(a) TECHNICAL PRESENTATIONS AND BRIEFINGS 


PI: N.S. GOEL 


P2: N.S. GOEL 


P3: N:S. GOEL 

D.E. STREBEL 

and 

R.L. THOMPSON 

P4: N.S. GOEL 
and 

R.L. THOMPSON 


"Progress in Scene Analysis Research: Use of 
Canopy Reflectance Models to Address Critical 
AgRISTARS Problems" Alov. 30, 1982, AgRISTARS 
Semi-annual Review, NASA-JSC, Houston. 

"The Use of Vegetative Canopy Reflectance Models 
in Support of AgRISTARS", Dec. 1, 1982, AgRISTARS 
Mini symposium, NASA-OSC, Houston. 

"On the Estimation of LAI and Other Agronomic Va- 
riables Using Canopy Reflectance Model, Especially 
The Suits Model". March 1, 1983, Quarterly Tech- 
nical Interchange Meeting, NASA-JSC, Houston. 

"On the Estimation of Agronomic Variables from Bi- 
directional Canopy Reflectance Data", April 29, 
1983, NASA-GSFC, Greenbelt. 


P5: N.S. GOEL 


P6: N.S. GOEL 

and 

R.L. THOMPSON 


P7; N.S. GOEL 


"Spectral Techniques for In-Si tu Measurement of LAI", 
May 1, 1983. Biological Productivity Experiment 
Review, HASA-JSC, Houston. 

"Estimation of Agronomic Variables Using Spectral 
Signiu.res", Sept. 12, 1983, Plenary Session. 

Inti. Colloq. on Spectral Signatures of Objects 
in Remote Sensing, Bordeaux, France. 

"Estimation of Agronomic Variables Using Canopy 
Reflectance Data", Oct. 19, 1983. Quarterly 
Technical Interchange on 'Assessing Key Vegetation 
Characteristics from Remote Sensing', NASA-JSC, 
Houston. 
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P8: N.S. GOEL 


P9; N.S. GOEL 
and 

R.L. THOMPSON 

PIO: N.S. GOEL 

and 

R.L. THOMPSON 
Pll: N.S. GOEL 


P12. N.S. GOEL 


P13: N.S. GOEL 


P14: N.S. GOEL 


"Estimation of LAI with Nadir View Spectral 
• Reflectance Data - Some Comments", October 
20, 1983, NASA-JSC, Houston. 

"CUPID Model and LARS Soybean Data", Oct. 20, 

1983. Workshop on 'Modeling of Crop Reflec- 
tance', NASA-JSC, Houston. 

"Optimal Illumination/Viewing Geometries for 
Remote Sensing of Agronomic Variables", Nov. 

7, 1983, NASA-GSFC, Greenbelt. 

"Agronomic Parameters from Bidirectional Canopy 
Reflectances", Jan. 10, 1984. Second Meeting 
of the NASA Fundamental Research Scene Radiatiore 
and Atmospheric Effects Characterization, Colorado 
State University, Fort Collins. 

"Agronomic Variables from Canopy Reflectance Data", 
March 9, 1984, Canadian Center for Remote Sensing, 
Ottawa. 

"Advanced Optical Techniques for LAI Estimation - In- 
version of Models", April 24, 1984. COVER Experiment 
Review on 'Characterization of Vegetation with Remote 
Sensing' NASA-JSC, Houston. 

"Optimal Illumination/Viewing Directions to Maximize 
the Accuracy of Estimation of LAI and LAD", April 
25, 1984. COVER Experiment Review on 'Characteri- 
zation of Vegetation with Remote Sensing', NASA-JSC, 
Houston. 


P15: N.S. GOEL "Estimation of LAI from Bidirectional Canopy Reflsc-,. 

K.E. HENDERSON tance", June 14, 1984. 10th Inti. Symp. on Machine 
and Processing of Remotely Sensed Data, Purdue University. 

D.E. PITTS 


(b) PUBLISHED REPORTS/PAPERS 


Rl: N.S. GOEL 

and 

D.E. STREBEL 


R2: N.S. GOEL 

and 

R.L. THOMPSON 


"Inversion of Vegetation Canopy Reflectance Models 
for Estimating Agronomic Variables I: Problem 
Definition and Initial Results Using Suits' Model". 
Remote Sensing of Environment, Vol . 13, Pages 487- 
507, 1983. 

"Estimation of Agronomic Variables Using Spectral 
Signatures". Proceedings of International Colloquim 
on Signatures of Remotely Sensed Objects. Bordeaux, 
France, Sept. 12-16, 1983. 
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R.L. THOMPSON 
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"Inversion of Vegetation Canopy Reflectance Models 
for Estimating Agronomic Variables II: Use of 
Angle Transfonns and Error Analysis as Illustrated 
by the Suits' Model". Remote Sensing of Environ- 
ment. Vol, 14, Pages 77-111, 1984. 

"Inversion of Vegetation Canopy Reflectance Models 
for Estimating Agronomic Variables III: Estimation 
Using only Canopy Reflectance Data as Illustrated 
by the Suits* Model". Remote Sensing of Environ- 
ment, Vol. 15, Pages 223-236, 1984. 

"Inversion of Vegetation Canopy Reflectance Models 
for Estimating Agronomic Variables IV: Total 
Inversion of the SAIL Model". Remote Sensing of 
Environment. Vol. 15, Pages 237-253, 1984. 

"Inversion of Vegetation Canopy Reflectance Models 
for Estimating Agronomic Variables V: Estimation 
of LAI and Average Leaf Angle Using Measured Canopy 
Reflectances". Remote Sensing of Environment. Vol, 16, 
Pages 69-85, 1984. 

"Estimation of Leaf Area Index from Bidirectional 
Spectral Reflectance Data by Inverting a Canopy 
Reflectance Model", Proceedings 1984 International 
Symposium on Machine Processing of Remotely Sensed 
Data. Purdue University, June 12-14, Pages 339-347. 

"Simple Beta Distribution Representation of Leaf 
Orientation in Vegetation Canopies", Agronomy Journal, 
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APPENDIX B 

FIRST PAGES OF THE REPORTS/ PAPERS WHICH HAVE BEEN PUBLISHED 

In this appendix, v/e have attached the first pages of the various reports/papers 
which have been published and which were listed in the preceding appendix. 
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Published Paper R1 

HiMon: sks’sisG oi HnvmosMEsr io - im - k ); udw ) -im 

Inversion of Vcgclalion Canopy llefleclancc 

Models for Eslimaling Agronomic Variables. I. 

l^roblcm Definilion and Initial Results Using the Suits Model 


NAREN’DRA S, COEL and DONALD E. STREREL* 

DiiwinHicin o/S|/j)o»w Srfpiff. S/Uff o/.Vnr VwL S'ew VoA‘ J3MI 


An ln)|>irUnt 1ml niitlvrl)' untnvvillt^-iU'd proli);in In rcinolc M'liilng h llii' Invctjlon of vcgclative tanop)* rcncclnnct 
inmIpU lo olilain iigtop}i)alcal iMtaintlcn, given nu'-vtnml iffli’cliinfei. Tlie prolilein fi licff lonnall)' dcfinwl and Hi 
ioltilion onlliwd. Nnnii'^'iil nonllnvar optlml^itlon lrtlinl«i't« we t»«l to Iniplcmenl iJili Invcrilon lo oblain llic leaf 
auM iiulcx ming Sulli' canopy tcflectancc model. Tlie leiiilti for a varfcly of Ciuci fndlealc llial lldi can lie done 
tnc'ceMlnll)’ iijliig fnftaietl lefk'CiiUicei at dlfletent cie'ii or lulmntlt anglei or a fombllWllon lliervof. Tlie dlier 
p.iriinielcrJ of Itie ino<iel imiil be known, iilibnngli le.iMmable inr,ijureinenl crrori can be tolerated witlimil cetfonily 
tlegradii'H llie acpiracy of tlic liiveriloii. Tlie applkallon of llie leeliiil(|tie to gionml Iviieil reinotc-ietiilng experinienli 
ix polcnUally melul, bm li limited lo ibe degtiv lo wldib the canopy lefleciancf model tan amiralely predict olnervcd 
lefleclancci. 


1. InlrocliicUon 

NNlieti elcclromagncllc |•.'l{liaUol) is iiid- 
dent on a vcgclalion canopy, il is scaUcred 
and reflected, and ils dircclion and spec* 
Iral composition arc altered in a complex 
manner by tlie vegetation. For llie pur- 
pose of ngncultiiral moivltoring by remote 
sensing, purl of Ibis altered and reflected 
radiation is intercepted by a salcllitc-bomc 
sensor. The success of Ibis monitoring 
melbod depends upon being able lo re- 
late reflectance measurements lo vcgc- 
lalion properties. One expects tbal Ibis 
can be done once the nature of tbe altera- 
tion in ibc radiation by the agropbysicnl 
iind environmental factors is specified. 

Over tlie last two decades tbcrc bau,ij 
been fairly intensive investigations, both 
experimental and lbcorclic.il, altcmpling 
to understand tbe rclatioasbip between 


*0n liMVe fioni Di'purtmvnl of 
Woffnnl Colli*)***, Sp.ir1.tnlmr};< SC 2^1301. 

» Stiriit v PitbliOiinj; Co » liR-.. 19H3 
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vcgclalton and environment variables and 
tbe .spectral reflectance. Tbese invcsllg.i- 
lions Itavc identified tbe key viiriablcs. 
On tbe tbcorctical side, one type of study 
has been the mo<lcliiig of crop canopy 
reflectance (see Goel, 1982; Smilb, 1982 
for reviews of various models). Here one 
defines or derives a function or an algo* 
rilbin wbicb yields tbe reflectance given 
llie variable,s .specifying tbe canopy. Tliesc 
variables generally include: optical prop- 
erlic.s of tbe vegetation components, c.g., 
wavelcngtb (X) dependent reflectance 
and traibsinitlancc of leaves, stalks, bends, 
etc.; physical parameters defining Ibc 
canopy geometry, e.g., density, angular 
iiicliuation, and dislribvillon of vegetation 
components; variables defining the soil, 
e.g., X-<lcpcndcnt reflectance of soil; vari- 
ables defining source of radiation (sun) 
and tbe properties of tbe detector, e.g., 
sun zcnilb angle, detector (ob.scrver) 
zcnilb and azimnlb (witli respect to tbe 
Min) angles. On Ibc experimental .side, an 
exlcn.sivc amount of data bas been col- 
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P ubllshed Pap er R2 


t 

ORlGeN/\l, P*'. ^ ^ 
OF POOR QUALh V 


ESTIILSTT ON O F ACHONOMIC VARIABLES 
USIKG SPECT.m SIGNATURES 


N.S. GOEL, R.L. THOMSON 
Dcpnrtncnt: of Systems Sclcnco 
State University of New York . 
Binskamton, New York, 13901 (U.S.A,) 


Sjy^hpjARY 

A technique for inverting a canopy reflectance model to estimate agronomic 
variables, using spectral signatures in the infrared region, is described. Its 
use is illustrated, for a soybean canopy, using a few representative models, due 
to Suits, Verhoef and Bunnik, and Norman. An analysis of the sensitivity of the 
estimated agronomic variables, like leaf area index and leaf angle distribution, 
to random errors in the canojjy reflectance data is presented. 


(Paper to be presented at the Int. Colloq. on Spectral Signature of Objects in 
Remote Sensing, Bordeaux, France. Sept. 12-16,, 1983.) 
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Inversion of Vcgclalion Canopy Reficclancc 

Models for Eslinialing Agronomic Variables. IT. Use of 

Angle Transforms and Error Analysis as Illuslralcd by Suits' Model 


I 

)! 

\\ 

•i 

ii 


NAHENDIU S. COEE, DONALD E. STHEBEL,* nnd mCHAHD L THOMI’SON b 

OiTWr/mtiif of S^tUvis Sciaicc, State of S\v> Yottc^ mnEjiarnton, AVu' Voffc I390! it 

ii 

Tilt u*clmiqtu? for liiviiilng a svgcIMlor) rAiiopy rcfkciiUKv itiodcl iwlkr (CocI and SUcl>cl» 1053) ii j| 

invciUgalcd hirtlKT. Tlic novel wnwpt of an ‘‘angle Imiijfonn** h Inlioduccd. Tills wnccpl allows Uie fonnalfon of j 

fnncUons of rcflccUnces si d^^’ '^nt vfew u*nlib and atfimitii angles, wldcli arc rlllit'r scnslllvc or fnscnsUfve to a I 

certain ngionomfc pirnmelcr. r r 0|K'r ('onibirmllon of ttiwc functions can allow delcnnlnalfon of all the Important 
agronomic and spectral paraine* ,r jom mc.wmcd canopy reflectance data. Tlie ietlmlcpif Is demonstrated ming Sulls‘ 

(1972) model for liomogenous canopies It Is sIjonsii tliat leaf area Index, leaf reflectance and transmittance, and overage 
leaf angle all can l>e determined from (lie tanopy reflectance at a set of selected Wew zenitli and azinmth angles^ A | 

sensitivity amdysls of the calndatin) values to tlic errors in tlie data Is also carried otit. Guidelines ore fonntdated for the i 

rmrnUr and types of observ'allmis required to obtain tlie values of n partladar canopy variable to within a given degree 1 

of accuracy for n given level of error In the measurement of canopy reflectance* I 


I* Tnlrodiiclion Improve the accuracy with wliidi tliosc jj 

paranujlevs arc dclcrnifncd wltcn mea- 

'llic cslfnmtlon of agronomic paranie- surcmcnl errors are present. i 

tors from rcflcdanco data Is an important We shall emphasize that an accurate 

practical problem. In the first p,apcr In esliinaticn of agronomic parameters from j 

this series (God and Strebcl, 1983} the meh^iircd canopy reflectance data has ! 

hereafter referred to as I), we defined this two key ingredients: 

problem formally and in some detail. 'Hie (1) A canopy reflectance model wbich ^ 

procedure for efiUinalion involves ins'cr- represents a canopy accurately, 

sion of a vegctallve canopy reflectance (2) A procedure for inverting such a 
model, 'rids inversion w.as Illustrated for model. 

two parameters, although it was noted In our initial studies, to focus our atten- 
that the use of more reflectance measure* lion on dcvcluplng techniques for invert- 

iiicnts would allow the determination of ing canopy reflectance models, we 

more parameters. Mere we extend the simulate the obsci'vcd canopy reflcc- i 

inclhod by showing how six reficclancc lances 'Dial is, we choose a ccrinin set of f 

data points, for different view directions, \’iilucs for the parameters occurring In a | 

can be used to obtain four parameters, model and use them to calculate the j 

and how additional data can be used to canopy reflectances. Tlie parameters ore 5; 

then ‘'forgotten,” and Oie calculated re* j 

*Oii t(Mve ftoMi ilir Di*p.taiiiftu of .si.itlicimittus, \Vo{. fluctniiccs arc taken as tlie observed ones. | 

lout Cotlvgc, Sp.iri.inlimR, sc 29301 . Tlicse ‘‘crror-frce” reflectances arc then ! 


-KlwUcf .Sdvmv I'liMIsliliig Co., )nc, lOS-l 
’)2 V,imlftl.iU Me, .S'fW York. NV 10017 
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Inversion of Vcgclalion Canopy Rcflcclancc Models 

for Eslimallng Agronomic Variables. III. Eslimalion Using Only 

Canopy Reflectance Data as Illustrated by the Suits Model 


NARENDnA S. COliL and RICHARD U THOMPSON 

Dt)»uf of Scfi’tiff, WflijoK SeUvoi of /riifi/iU'rrinfc Sloic 1‘iift tvif/v of .Vnr Vort, Bnpg/iCHifort, ,VV 1390J 


Tlie tcclinkjiic lor villiiiiilitig agronomic and ipcclriil lutiinivicrt for a M-gcialloti c.uiop>' from llic canopy tcncclaiicc 
(CU) ctaln, (Jrscrilrod earlier (Coel el nl„ I0S4), liiit been liiiproveil Tliwe ImprovemeiiU are Iwofold: Iltri, one wn 
now, in prindplc, cJliinale varlom p.iraincler> ming only CR data, In lire Intrared band, lor a »ct ol »olar/>1ew angles! 
secotrd, tlie mclbod is now compulatlonally nineb more ellldent. Tliese improvenieiUs are lllnstroli'd sia Suits' mrnlel. 
An aiiidysts of lire sensitivity of the calcidaled agronomic mid spectral luraniclers to ebanges In ibe CR is also i.nried 
out. lids analysis suggests llinl, In general, for cs|ieclcd levels of errors In lire measurement of CRs and tlie atviiracy 
Mitb which the Suits model Is lihely to represent CR, one Is unlihely to Ik> nhle to estimate agionoinie parameters like 
leal area Index (LAI) oiid average Icid luigle (ALA) using only measured CR data. Such a determination will likely 
reipiire aiieill.iry data on the reflectance and tran>iiiiliance ol vegrtallon elements and on the soil reflectance. 


1, Iitlroduclion 

One of Uic most desirable goals of tlie 
research on the remote sensing of vegeta- 
tion canopies is to be able to estimate key 
agronomic parameters like leaf area index 
(LAI) and leaf angle distribution (LAD), 
using only canopy reflectance (CR) mea- 
surements. LAI is perhaps the most im- 
portant variable for dclcnnining growth 
and yield. 'Jlioiigh LAD is not of the 
.same importance as LAI, it is e possible 
indicator of plant stress level. With the 
goal of e.stiinaling agronomic parameters, 
many systematic theoretical and experi- 
mental invcsllgntlons have been carried 
out in the last quiirtcr of a cenltuy. These 
lnve.stigations have led to the clarification 
and quantification of the relationship 
between these and other p.il ,tmelers de- 
scribing the spectral properties of vt-gcla- 
tion elements and soil, incident sokir flux, 
and viewing direel ions iind the canopy 

' SifoMtv PiililiOihit; Co., Inc., J9S*1 

52 N'imik*fblll MVu New York, NY HK)17 


icfleclancc in various Wiavclcngth bands 
(Ross, 1981; Goel, 1932). 

Our Invcsllgalions have focused on de- 
veloping an approach for obtaining agro- 
nomic paramclcrs by inverting a CR 
model which rcpre.scnls the relationship 
between canopy paramclcrs and CR. For 
this purpose wc have Initially chosen the 
Suits model (Suits, 1972) for a homoge- 
nous canopy. In the first paper in tills 
aeries (Goel and Slrebcl, 1983; hereafter 
referred its I), wc defined this problem 
formally. We showed that one could, in 
principle, estimate LAI and average le.af 
angle (ALA) using CR data, in the in- 
frared iiaiid, for seven'll view angles, pro- 
vided one knows the otlicr auxiliary 
parameters like canopy component (c.g., 
leaf) rcflcclancc p, and transmittance t, 
.soil rcfleelancc p,, and fraction of dif- 
fused incident solar radiation, SKVL, In 
the .second paper in this scries (Goel et 
al, 1984; hereafter referred as II), wc 
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iiFMon: s/,^'^/yc on:s\’mw\u:sr i5:iM7-a53(jos4) 

Inversion of A^egcUUion Canopy Rcfleolancc Models for 
EsUnialing Agronomic Variables. XV, Total Inversion of the 
SAIL Model 


NAHENDIU S. GOEL mid RlCJlABD L. THOMPSON 

DiV>ift»ii m of S|/iln)u Sfimcf, M’rtficm Sf/ioo/ of S'lCfm’rWtic, SJufr ('iificnffu of S'ni' YoiK Binef>otnton, AT J3007 


Tlic wllmabllfly of stt llic evinop)' piifamcli-rj (or a vcgclallon c.iitopy uiliig only canopy icnccIniHc (CR) data iiml llic 
SAIL model is Invcsllgalfd, usIiib Icclmiqncs dcstriltcd earlier (Goel and Tliompson, lOSia). Il Is shown dial In 
principle siicli an esllniailon is i«sslWc, i c„ llic SAIL morlcl is inallieniallcaHy lolaliy Invertible. An 'iunlysls ol llic 
sensllisily ol llie calculaled values lo cliaiiges In the CR data is presented. Tliii analysis siigficsts »hnt, given the 
expecled atcnracy o( CR ineasutenienls and the acetitaey ol the .SAIL model in tepreseiillns CR In the inirared region, 
Ihc agronomic paranielers,leai nrea Indea, and leal angle dislrihiitioii, can be csllmaled fairly accairatcly using Hiicillary 
data on spectnd parameters. 


1. Inlroduclion 

Re I Hole sensing of vegclnllon relies on 
the speclnl aignnlnres of the vegetation, 
in parlictiliir, the signatures in the visible 
and infrared regions. During the last 25 
years or so, many models (Goel, 1982j 
Ross, 1981) at different levels of complex- 
ity have been developed. They nltcmpl to 
provide a realist ic relationship between 
the important agronomic and spcelra! 
parameters of Uic vegetation canopy and 
its canopy reflectance (CR) signatures. 
The agronomic parameters arc leaf nrea 
index (L.AI) and leaf angle dustribitlion 
(LAD). The spectral parameters include 
leaf reflectance p and transmittance t. A 
desirable application of these models is 
the estimation of these parameters from 
the measured rcfleclanoes. 

In this .scries of papers (Goel and 
Strebcl, 1983; Goel et al, 1984; Goel and 
'niomp.son, 1984a) \vc Iwive been investi- 
gating the possibility of such an estima- 
tion, using CR data in the infrared region 
for a set of solar/vicw directions'. In the 

* Rlsi’VltT Stii'iuv Piililidiliig Co., Im*.. iOS4 
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preceding paper of this series (Goel and 
'Hiompson, 1984a; henceforth referred to 
as 111), we have shown tliat such an 
estimation is possible, in principle, for the 
Suits model (Suits, 1972) for a homoge- 
nous canopy. The parameters occurring 
in this model are the leaf reflectance and 
transmittance, leaf area index and aver- 
age leaf inclination angle Oj, soil reflec- 
iance p,, and the fraction of diffused inci- 
dent solar radiation (SKYL). Specifically, 
we showed that all of these parameters 
can be estimated using only CR data and 
the Suits model, provided that the Suits 
model is an accurate representation of 
canopy reflectance and that the lencc- 
lance measurements are accurate. In other 
words, the Suits model is mathematically 
totally invertible. We also analyzed the 
.sensitivity of tlie parameters estimated by 
model inversion to changes in the CR 
data. 

In this paper we apply the techniques 
developed in 111 to a more complex model 
—the so-called SAIL (scattering by arbi- 
trarily inclined leaves) model due to 

00W-i2.57/H/$3.(K) 


0 RI 6 INAI. l:Q 
OR POOR QUALIW 


GO 


-Ct- 

Published Paper R6 

HIMOn i>KSSISC 01- L'.W'inOXU/iSX IG^GO-fi-J (JfJS J) 

Im’crsion of Vegclnlion Canopy Kuflcclancc Models for Eslimaling 
Agronomic Variables. V, Eslimalion of Leaf Area Index and 
Average Leaf Angle Using Measured Canopy Refleclanccs 


.NARENDRA S. GOEL AND RICHARD L. THONU’SON 

nkiiurfim'iif D/Syirciiii Sffnifc, U’d(»wi .Sr/ioo/ o/'CiifiJiU'crfitc, Stair t'ltlim/fKo/.Vni’ Vorfc, M' 1390J 


Till' li'ilmiijHi' (Imtilicil imfIIit (Gwl mid Tliomp>on, 19Sa)) lor cMlmnlliiniiBinnmiilc pirmnclcrj from lildlrtclioiwl 
rtop ri'ni'tlmuf ftiila iJ (ippHi'd to o fully iwcri'd wylm.in liiliopy, lining ilntn mcaiiiri-d In llic floW. Tliii U'clinltiue 
oinployi llic liivi'tilon of a c.inopy ii'fli'i'lima’ model it f> >lioin> tliiil luiiig llie SAIL moilel one c.ui eilfmote leaf area 
Imle* (LAI) aa well as as er.ige Iwf angle (ALA) ipiite well, proeliled llial Hie oilier iiuiopy paraiiu’ie>« (leaf reflectmiec 
iiml l!‘.iii4tnlllaiicf, anil tonectame, and fr'ieUoii of diffmed il>yHg)it) are litiouii. Some siiggeslloiu are made for 
Impioeing llie S AIL iiioilel Tlili dionid Improve the aeenraey of edimallon of nol only LAI and ALA tmt sliould also 
allow the eMImatloii of |)ie eompUie leaf angle dlilrlluitinn. 


1. Inlroduclion 

It is now vcjy well estiTbllsbcd Ihal the 
R'llcdiinccs, in (lie visible nnd infrared 
regions, from a vegclnlion canopy arc 
slrongly corrcl.ated lo (he agronomic and 
spcclral parameters of (lie canopy. The 
agronomic pnramelcrs include leaf area 
index (LAI) and leaf angle dislribution 
(LAD), Tlie spcclral paramclers Include 
leaf reflectance p and Iraiibniiltancc t. 
Tliis eonelalion has been quanllfied by 
many canopy reflcclancc (CR) models 
(see Goel, 1981; Smilb, 1983; Ross, 1981 
for reviews of some of these models). It 
will be vcr>' desirable to exploit this 
coirelnlion to estimate llie .agronomic 
pammelcrs from the measured canopy rc- 
fledimcc data. 

In this sciles of piipcrs (Cod and 
Slrcbel, 1983; God cl al., 1084; Cod and 
Thompson, 1 984a, b) we have been In- 
vest igaling the possibility of such an 
eslimalion, using CR data in the infrared 
region for a set of solar and view direc- 
tions. Wc have shown llijtl such an esli- 
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maiion is, in principle, possible al least 
for a homogeneous vegetation canopy. 
The procedure involves the inversion of a 
canopy reflectance model. In the prccccd- 
ing paper of this series (God and Tliomp- 
son, 1984b) and in another paper (Cod 
and Tliompson, )9S3) we showed that 
one can determine all the canopy param- 
clcrs occurring in the SAIL model 
(Verhoef and Bunnik, 1981) and in the 
CURIO model (Norman, 1979) using only 
CR (lala in the infrared region. The 
paramders occurring in these models are 
the leaf reflectance and lransmiUance,leaf 
area index, and leaf inclination angle dis- 
Irihulion (LAD), soil reflectance p^, and 
the fraction of diffused incident solar 
radiation, SKYL, Thus all of these param- 
clcr.s can be estimated using only CR 
data and Ihc.se models, provided llial these 
models aecuralely represent measured 
cimopy rcfledanccs. 

To focus our initial allention on dcvel- 
opitig Icchniqucs for e.s(imating canopy 
parameters, wc have so far simulated tlie 
ob.sem'd canopy reflectances in our stud- 
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ABSTHACi' 

A technique for estlmu t Inc tkc leaf 
area Index from bidirectional canopy re- 
flectance (CR) data. In the Infrared re- 
gion, in band 4 of a Thematic 

Mapper (TM), is described. It involves 
inversion of a CR model which accurately 
represents the reflectance from the canopy, 
A method for remotely collect Inc this CR 
data using an aircraft based TM is de- 
scribed. The bidirectional CR*s, for a 
black spruce (picea naricna) canopy, for 7 
solar/view directions, as measured using 
this technique, arc given, A very pre- 
liminary analysis of the data from a point 
of view of ostlrnatinc LAI by Inversion of 
r. CR model Is given. This analysis 
sucCGsts that for an acceptably accurate 
estimation of LAI, one will require bi- 
directional CR‘s for many more than 7 
solar/vlew directions. 

I. INTRODUCTION 

It is now very well established that 
the reflectance, in the visible and in- 
frared regions, from a vecetatlon (crop, 
Srassland, and forest) canopy is strongly 
correlated to the acronomlc, archi tectuial , 
ar.d spectral parameters of the canopy. 
Remote sensing of vegetation relies on 
this correlation. The agronomic parameto*s 
Include the densities and orientations of 
vegetation components like leaves, stems, 
branches, and bark. The architectural 
parameters Include spatial d is t r I bu t ic»ns , 
both in horizontal and vertical direc- 
tions, of vegetation components. The 
spectral parameters include reflectances 
and t ransmi t tances of these components. 

In addition, the canopy reflectance (CR) 
depends upon the ground (soli, moss, 
disintegrated vegetation etc.) reflectance 
and the relative fraction, SKYL, of the 
diffused incident solar radiation. 

During the last quarter of a century, 
many models (sec Ross, 1981; Gocl, 1982; 


and Smith, 1983 for reviews of these 
models) at different levels of complexity, 
have been proposed to provide a realistic 
relationship between important canopy 
parameters and canopy reflectance. With 
these CR models one can, in principle, cal- 
culate the canopy reflectance as a function 
of the illumination and view directions 
(bidirectional CR) , using the measured 
values of the canopy parameters, ground 
reflectance, and Incident diffused radia- 
tion. A desirable application of these 
models is the ostlmatlpn of important 
agronomic parameters like leaf area index 
(LAI) from the measured bidirectional CR’s 
by carrying out such calculations in 
•‘reverse**, i.e., by inverting a CR model. 

One of the authors and his collabora- 
tors (Goel and Strobel, 1983; Cool, 

Strebel and Thompson , 1984; Goel and 
Thompson, 1984a, b,c) have investigated the 
possibility of such an estimation, using 
CR data in the infrared region, for a set 
of solar/vicw directions. We have shown 
that such an estimation is, in principle, 
possible, at least for a homogenous 
canopy. Specifically, we have shown that 
at least for three CR models - the Suits 
model (Suits, 1972), the SAIL model 
(Verhoef and Bunnik, 1981), and the CUPID 
model (Norman, 1979), there is a one to 
one relationship betw’een the canopy para- 
meters and the bidirectional CR. That is, 
all the canopy parameters can be estimated 
uniquely using only CR data, provided of 
course these models accurately represent 
the measured CR's. V/e also applied tL..s 
canopy reflectance model inversion tech- 
nique to field measured CR’s for a set of 
view directions. V/e have shown (Goel and 
Thompson. 19S4c) that one can estimate 
quite accurately the leaf area index (LAI) . 
as well as the average leaf inclination 
angle (ALA) for a homogenous fully covered 
soybean canopy^ using CR data for about 50 
view directions, provided that the other 
canopy parameters like leaf reflectance 
and transmittance, soil reflectance and 
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Simple Be(a Disliibution RepresentiUlon of Ij?af OrionlnUon in Vegetalion Canopies' 

Nnrendra S. GocJ nnd Donald E* StrebeP 


ABSTRACT 

It is slifmn tltaj 0 tviO'parLmcItr beta distribution rt-jircsenls ^cr)‘ 
>u*ll the leaf nniilc- distribution for o \nr|ely of )C);cl;ttl(in c.inopics, 
Compailsun to lhturciir*il distributions as noil us field data sii^^ests 
that this diviribulifin be runsidcfcd ns « candidnic fur a “unlursal” 
leaf un;;le Histriholinn, The f urRincters of the dislrUmtlon c;m then 
bo used to Idcntlf) and classify Initrinedinte canopy <>pos mid pos* 
siblj to qiumlfy ciimip) chun};es mused b) cn^liunmontul stress* 

4dtIifhncJ hdt'jc uardsf Uaf anylc distribution, Crop clussUim* 
tton, Remote sinslog. 


' Phe arrangement ofleavcs in a plant canopy influ- 
A cnees the Interaction ofclectromagnctic radiation 
with the plants. In particular, photosynthesis and yield 
arc determined by the light distribution within the can- 
opy and hence by the a?iimJth, inclination, and spatial 
distributions orieavcs. These distributions also deter- 
mine the ladlation reflected from the canopy, thus in- 
fluencing the spectral signature used for the purpose 
of remote sensing (Smith, 1983). It is iheicforc im- 
portant to liavc a good way to describe the leaf dis- 
tribution. 

The distribution of leaf inclination and azimuth .an- 
gles vary fiom crop to crop and may also be dependent 
on lire growth stage of the crop (Ross, 1981) or e\cn 
the time of the d.iy (Kimes and Wrehner, 1983), One 
may CApeci changes ;d so to be correlaicd with stress 
.'ind hence changes in the.se distributions could be used 
to ()ii:inlify ihc .suess. Tor species uiih no preferred 
azimuth direction, dc Wii (1965) and others (Russ, 
1981; .Smith, 1982) ha\e developed scwtal special dis- 
tribuiions to ch.:ntcierizc the leaf angle disirilniiion 
(LAD) uhich distinguish m.ijor kinds of pl.mi c;mo- 
pies. These tspesarc: (l)planophilc— hori/onial leases 
most freijueni (2)t rrclophilc— s'crtical leaves most fre- 
<iOem, (3) plagiophilc -oblique leaves most frequent, 
(4) eviremoj l.ilc -oblique leas'es le.sst frecjiient, (5) 
unifoim -vaine proporiion of le.:ve5 at ans ang.le, and 
(6) splii i ii ;d -hvif angle frequency same as for suif.tcc 
eleiiu nts i>f a jjiherc. This case b) case apj)ro.'»ch, 
though is sumewh.nt limited in analszing inletmediaic 
canops tspes or the i.’inre of s ariation within a c.mopy 
i>pe. 

In this intie sso po',” >se that the sselld.nossn two 
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